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Abstract 

Concepts envisioned for the future of civil air transport consist of unconventional 
propulsion systems in the close proximity to the structure or embedded in the airframe. While 
such integrated systems are intended to shield noise from the community, they also introduce 
new sources of sound. Sound generation due to interaction of a jet flow past a nearby solid 
surface is investigated here using the generalized acoustic analogy theory. The analysis applies 
to the boundary layer noise generated at and near a wall, and excludes the scattered noise 
component that is produced at the leading or the trailing edge. While compressibility effects 
are relatively unimportant at very low Mach numbers, frictional heat generation and thermal 
gradient normal to the surface could play important roles in generation and propagation of 
sound in high speed jets of practical interest. A general expression is given for the spectral 
density of the far field sound as governed by the variable density Pridmore-Brown equation. 
The propagation Green’s function is solved numerically for a high aspect-ratio rectangular jet 
starting with the boundary conditions on the surface and subject to specified mean velocity 
and temperature profiles between the surface and the observer. It is shown the magnitude of 
the Green’s function decreases with increasing source frequency and/or jet temperature. The 
phase remains constant for a rigid surface, but varies with source location when subject to an 
impedance type boundary condition. The Green’s function in the absence of the surface, and 
flight effects are also investigated. 


1. Introduction 

| OMMERCIAL aviation in the U.S. is projected to increase by 70% in the period 2010 - 2030, resulting in 
^■^more exposure to air traffic noise in the communities surrounding airports. Under the Fundamental 
Aeronautics Program at NASA, the Fixed Wing (FW) Project is tasked with exploring concepts and technologies 
that could significantly improve energy efficiency and environmental compatibility of the fixed wing subsonic 
transport aircraft [1]. FW research focuses on vehicles that are three generations (N+3) beyond current state of the 
art (N) aircraft, with noise reduction goals of -52dB by 2025-30 time frame (referenced to 737-800 aircraft with 
CFM56-7B engines). In 2008, NASA instigated a call for conceptual design of future commercial transport in order 
to meet specific goals related to noise, air pollution, and fuel consumption. Early studies indicated that both 
propulsion system and engine placement have to be modified in order to achieve these new objectives. In particular, 
placement of the aircraft engine for the purpose of reduced community noise poses new challenges in structural 
design, material selection and computation. Configurations such as over the wing engine mount, distributed 
propulsion, Hybrid Wing Body (HWB) concept, and/or high aspect ratio rectangular exhaust geometry with 
extended beveled surfaces shield noise from reaching the ground. These concepts also produce new noise sources 
due to scrubbing of the flow past structure, skin vibration, boundary layer (BE) noise, and edge noise due to the 
scattering of sound from sharp edges. 
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In 2010, a series of tests were initiated at NASA Glenn Research Center in order to study the propulsion/airframe 
integration. These experiments, under the umbrella of “Jet-Surface Interaction Tests (JSIT),” were intended to guide 
analytical studies towards development of prediction model. The initial tests consisted of placement of a flat plate in 
the proximity of a round jet. Details of the experimental setup and test cofigurations are provided by C. Brown [2], 
and G. Podboy [3]. Measurements were gathered on both sides of the plate, commonly referred to as the reflected 
and shielded sides. In addition to the far-field microphone noise measurements [2], phased array source localization 
experiments [3] were also carried out in order to shed light on the noise source location and its frequency content 
under different exhaust conditions. 



Figure 1. Sample JSIT noise measurement , C. Brown [Ref. 2]. 


An example of JSIT noise data that 
shows power spectral density per 
Strouhal frequency in a 2-in diameter (Z) 
=5. 08cm) Mach 0.97 unheated round jet, 
as reported by Brown [2], is shown in 
figure 1. Here a solid surface extending 
12-diameters from jet exit in the 
downstream direction was positioned at 
DU from the jet centerline. Far-field 
measurements are shown at 90° and at 
100Z) for a jet in isolation (no solid 
surface), as well as shielded and reflected 
sides of the jet in the presence of the 
surface. Examples of source localization 
data for the above jet are reported by 
Podboy [3] (not shown here). These 
measurements present beam-form maps 
of the source location for an isolated jet 
as well as when surfaces of varying 
length are placed between jet and the 
phased array. 
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Figure 2. Jet noise spectral measurements at 90° in a Mach 
0.97 rectangular exhaust, J. Bridges [Ref. 4]. 


Similar measurements were reported by 
Bridges [4] using rectangular jets of 
varying aspect ratio (AR) in the 
proximity of a flat surface. Sample 
spectra are shown in figure 2 for a Mach 
0.97 unheated rectangular jet (AR=8), 
with dimensions of (5.3566 x 0.6696-in), 
and the major axis parallel to the plate at 
standoff distance of h = 0, and plate 
length of X TE = 12.0-in. Measurements 
are shown at a polar angle of 90°, i.e., 
directly above or below the surface, as 
well as for a jet in isolation. An 
extensive database has been generated as 
a function of the jet operating conditions, 
length and position of the solid surface, 
and observer angles. 
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It is commonly accepted that the trailing edge noise is responsible for the enhanced low-frequency and peak noise 
amplitude measured on both reflected and shielded sides of the plate as seen in figures 1 and 2. A distinctive feature 
of the trailing edge noise component is its cardioid directivity shape and a dependence on jet velocity to the fifth 
power [5,6]. Recent measurements by Lawrence et al. [7] indicate that the velocity power factor could be as large as 
6 at longer plate length L that yield a wetted trailing edge. FFowcs Williams and Hall [5] describe a solution to the 
edge noise from a rigid plane in the presence of a turbulent flow using Lighthill’s acoustic analogy. They consider a 
very low Mach number flow, and take advantage of an appropriate Green’s function applicable to a point source in 
the vicinity of a half-plane [8] . Howe [9] reviews various theoretical approaches in understanding the trailing-edge 
noise, and summarizes diverse theories into a single model that exhibits the dependence of the acoustic field on the 
turbulent fluctuations near the edge as well as the significance of the Kutta condition at this point. Goldstein et al. 
[10, 11] include the effect of the mean velocity gradient on the trailing edge noise by imposing a convecting vortex 
source, or gust, of the general form / U (x 3 ) — t,X 2 ,x 3 ) on a transversely sheared mean flow. 

Our interest here is the turbulence-generated noise due to the scrubbing of the jet flow past an adjacent flat surface. 
At very low Mach number, Howe [12] describes the far- field acoustic domain in terms of the wall-pressure wave 
number-frequency spectrum. Various models have been suggested for wall-pressure spectrum [13] that are used in 
conjunction with a 2D or 3D free-space Green’s function to express the far-field spectral density. These 
approximations ignore the mean flow refraction, compressibility effects, and sound produced outside the BL at 
distances that are large relative to the acoustic wavelength. At jet velocities of practical interest typical of a jet 
exhaust, frictional heat generation and thermal gradients in the BL could play important roles in generation and 
propagation of sound. Additionally, skin vibration could influence the surface boundary conditions and the 
subsequent radiated sound. 

In sections 2 and 3 we derive an expression for the propagation Green’s function (GF) applicable to the scrubbing 
noise on a flat surface. The GF needs to be obtained numerically when the mean flow profiles are non-uniform. An 
expression for the far-field pressure is given in section 4. Ideally a Reynolds-Averaged Navier-Stoke (RANS) 
solution for the jet in the proximity of the surface provides the mean flow and turbulence information required in 
evaluating both the source strength and the GF. A 2D approximation of the GF (when the flow is infinitely long in 
transverse X 2 direction) is discussed in section 5. Sample GF computations using analytical profiles for the mean 
flow and a summary are provided in sections 6 and 7. 


2. Formulation of the Scrubbing Noise Problem 

Consider turbulence-generated noise due to the scrubbing of a jet flow with an adjacent flat plane (Fig. 3). 
The governing acoustic perturbations may be obtained by writing the Navier- Stokes (NS) equations as a set of mean 
flow equations (referred to as a non-radiating base flow), plus a set of linear equations for the fluctuating 
components of the motion. By selecting a set of five appropriately defined perturbation variables, four of them non- 
linear, Goldstein [14] shows that the left-hand side of the acoustic equations resemble those obtained by linearizing 
the convective form of Euler equations about a similar base flow. Further, by assuming the mean flow as locally 
parallel, the factors on the left-hand side of momentum and energy equations that explicitly depend on the viscous 
stresses reduce to higher order terms and are neglected. These equations would then convert to inhomogeneous 
Rayleigh equations where the viscosity effect only appears as a non-linear source term. 

In formulating the scrubbing noise problem, the mean flow is represented as a two-dimensional sheared flow. We 
consider the mean static pressure p within the BL to be a constant. At very low Mach numbers the compressibility 

effects may be negligible, however at Mach numbers typical of a jet exhaust, frictional heat generation and thermal 
gradients in the BL could play important roles in generation and propagation of sound. The RANS solution (which 
is now considered as the base flow) provides such effects via the action of viscosity near the surface. Hence by 
maintaining the compressibility effect, noise generation due to the scrubbing of a jet past a nearby surface is, in 
general, described by full acoustic analogy as considered in Appendix- A, and by Pridmore-Brown equation within a 
locally parallel mean flow approximation [15] 
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The pressure variable Tl' consists of additional (higher order) terms that reduce to zero in the far field acoustic 
domain, and the source term F(x,^) on the right hand side of (1) is defined according to the generalized acoustic 
analogy [14]. 
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stress tensor is the difference between the fluctuating and Favre-averaged Reynolds stresses 


e^-ipv'v'-pv'v'), 


and /r denotes the fluctuations in the moving-frame stagnation enthalpy (see the nomenclature). 


(3a) 


(3b) 


(3c) 



Figure 3. Jet scrubbing on a nearby surface. 


The mean velocity U and density p are functions of the normal coordinate (jc 3 > 0), and c(x 3 ) is the mass- 
averaged speed of sound 

U = U(x 3 ), c 2 (x 3 ) = yp / p(x 3 ) . (4) 

Note that, in practice, both the mean flow and turbulence depend on the transverse direction x 2 as the source strength 
eventually diminishes at some distance ±L 2 away from the jet centerline. In the stream- wise direction, the 
scrubbing effect may be limited to a distance 0 < X l < L x , or the jet may persist beyond the plate edge at L\ where 
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interaction of the turbulence with the trailing edge produces edge noise. Here, in formulating the scrubbing noise 
problem, we limit the source to a finite x 2 span by considering the mean flow as locally parallel in both x\ and x 2 
directions - thus expressing the problem in three dimensions. This assumption, although more suited to jets with a 
high aspect-ratio rectangular exhaust, permits source- volume integration within the wetted span in the x 2 direction. 
A round jet may be transformed into a rectangular strip by using an appropriate conformal mapping transformation. 
When all changes in span-wise x 2 direction are neglected, the sound intensity should be expressed as per unit 
thickness. The special 2D case will be discussed in section 5. 

We seek a solution to equation (1) of the form 


n\x,t) = ^G(x,t\y,T)T(y,T)dTdy , (5) 

y T 

where G denotes the GF 

LG(x,t\y,T)-8{x-y)8{t-T) . ( 6 ) 

In a bounded medium, the above form of the solution implies that the GF satisfies certain boundary conditions on 
the surface. The source- volume integration in (5) is carried out over noise-generating region y in the vicinity of the 
scrubbing surface and at the emission time T . The solution described here excludes the trailing edge noise that 
travels to the opposite side of the surface as alluded to earlier. This latter noise component requires a Wiener-Hopf 
type solution. 

We define a Fourier Transform (FT) pair with respect to (x x ,X 2 ) coordinates as well as time t. The corresponding 
wave number vector is denoted as k t — (k x ,k 2 ) 

G(k t ,x 3 \y 3 ,a>)= ^ J j j G(x4\h^) e ~ ,Wx '~ yi ^ X2 ~ y2ya °~ T)] d(x-y l )d(x 2 -y 2 )d(t-z) , 
G(x,t\y, t)= J J J G(k, ,x 3 \ y 3 ,co)e i[ k] ’V/A’, J/: 2 . 


Throughout, a circumflex is used to denote a FT variable as defined above. 


The transform of equation (6) is 


d 2 G 

dxl 


- + 


(c 2 )' 2 kJJ' | dG 

c 2 -co + kJJ J dx 3 


- + 


(-co + k.Uf 


G - 


Hx 3 -y 3 ) 

(2 nf Si-co + kJJ) 


( 8 ) 
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where prime notation, as used on mean variables C and U, signifies differentiation with respect to normal distance 
X 3 . The equation possesses a critical point atx 3 where factor —CO + kJJ (x 3 ) vanishes. To avoid this singularity 

we require k x < CO I U e where U e denotes some peak mean velocity in the BL. This upper limit for k\ is usually 
referred to as the non-convective domain for the wave number. As discussed in [16, 17], in the convective domain 
k x > CO IU e it becomes necessary to provide a loss mechanism, such as turbulent viscosity, in order to avoid the 

singularity. The dominant value of k x contributing to the far- field noise will be determined in section 3. It will be 
shown that when U e < C ^ the upper limit condition stated for k x will be satisfied at all observer location. 

For convenience, the second-order linear differential equation (8) is rearranged as 
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(9) 


d 2 G 

dx] 


+ f(k t , x 3 ,co)G = 8{x 3 - y 3 ) . 


Since equation (9) is self-adjoint, the newly defined GF, G , is symmetric in (x 3 ,_y 3 ) 
G(k t ,x 3 1 y 3 ,co) = G(k n y 3 lx 3 ,ft)) , 
and is related to the GF of interest G as 


G(k t ,x 3 \y 3 ,co) = 


1 -co + k l U{x 3 ) 


{Inf c(y 3 )c(x 3 ) {- (0 + k x U{y 3 )) 


G(k t ,x 3 \y 3 ,co ) , 


where 


/\2 


f(k t ,x 3 ,co) = x 2 - C 2 2{k ' U) +. , 

; 3 {-co + k.Uf -co + k x U c 2 -co + k.U 

X 2 {k,,x 3 ,co) = - — — kf-k 2 . 


kJJ" { c 2 )' k x U' 
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i a f(c 2 r) i^- 2 ^ 2 


2 dx. 


(c z y 


c y 


+- 


(c 2 y 


V c 
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(id 


and (c 2 )' / c 2 = {dp / dx 3 )/ p . 

It should be noted that unlike the present problem, the full propagation equations (i.e. not locally parallel) might not 
reduce to a self-adjoint form. Such is the case when the governing equations are written for a spreading flow [18]. 
In that case, working with the adjoint GF will entail additional integrals that should be carried out on the surfaces. A 
general solution for the five field variables using a vector GF formulation is provided in Appendix A. The second 
integral in equation (A 16) in the appendix describes the contribution to the acoustic field due to the nearby surfaces. 
In the formulation described herein, the GF satisfies conditions that make the surface integral equal to zero. 

Far away from the noise generating region, as X 3 — > 00 the mean velocity and the sound speed approach their 
respective ambient values U ^ and C ^ , and equations (9) and (11) show that 


= x 3 ->°o 

where 

X 2 ={-k+LM ) 2 -k 2 -k 2 \ K=(0/c , M =U Ic . 

/t 00 \ O foo/ 1 Z ' O 00" 00 00 00 


(12a) 


(12b) 


The Mach number M ^ takes account of the flight effect when surface is in motion relative to the ambient. The far- 
field solution to equation (12a) is 

G(k t ,x 3 \y 3 ,co) = b e~ lXo ° x \ x 3 — > °° (13) 


where complex number b has dimension of length, and depends on k t and (0 . An outgoing propagating wave ( 
X 3 > 0 ) requires the negative root of xlo when > 0 . When xto < 0 the pressure decays exponentially 
normal to the surface (evanescent waves), and the branch-cut is chosen such that Xoo — — Z*0 I) 12 • Aerodynamic 
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noise scattered off the edge of a half-plane is sometimes described as the interaction of evanescent waves (rather 
than propagating acoustic waves) with the edge [19, 20]. These disturbances are dominated by spectral components 
that possess a subsonic phase velocity in a constant x 3 plane. 


Since we are interested in the scrubbing noise that reaches a distant observer, the radiation field would be dominated 
by wave numbers, in constant x 3 plane, that possess a supersonic phase velocity 


CO 


M 


-K + 




> c„ 


(14) 


A no slip boundary condition on the mean flow (U = 0,X 3 — > 0) combined with the x 3 component of the 
linearized momentum equation within the turbulent boundary layer (TBL) shows that 


d(Pv j) t a / = Q 

dt 3x, 


(15) 


Upon replacing p with the mean fluid density p o near the surface, we apply a FT to the latter equation to obtain 


O dp' 

-i(0p o v 3 + — = 0 , 

ox. 


X 3 “ 0 


(16) 


Assuming that Fourier components of the normal velocity and pressure on the surface relate through the surface 
impedance Z(co) as 


we find 


p = Z v 3 , x 3 = 0 


= 0, x, = 0 . 

dx 3 Z 3 


Since the FT of a convolution is the product of the Fourier transforms, equation (5) implies 


(17) 


(18) 


7t'(k t ,x 3 ,(0 ) = (Inf ^G(k t ,x 3 1 y 3 ,co)T(k t ,y 3 ,0))dy 3 , 

y? 


(19) 


and upon using (19) in (18) we find 


dx 3 Z 


X 3=0 ■ 


( 20 ) 


Our interest is in the boundary condition applicable to G rather than G . Substituting equation (10) in (20) shows 
that 


y/G-0, x 3 -0 , 

OX 3 


(21a) 
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yf(k l ,(o,Z) = 


(21b) 


i K c 2 

0 _ °° 

~Y c 2 


c' k ^ 

+ ^ + ^U'( 0 ) , 

C„ CO ) 


Z(co) = 


Zico) 
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where C o is the mean sound speed in the fluid at the surface interface, and C o is its derivative at x 3 = 0. An 
alternative way of deriving the surface boundary condition (21a) is described in Appendix B. 

The problem at hand is thus reduced to solving equation (9) subject to the two boundary conditions (14) and (21). 


The GF in equation (9) may be expressed in terms of two linearly independent solutions Vj(k t ,X 3 ,(Q) to the 
homogeneous equation 


d 2 V j /dx 2 3 +fV j = 0, j - 1,2 


( 22 ) 


such that Vi and V 2 satisfy the homogeneous boundary conditions at 0 and oo, respectively. For brevity we omit k t 
and CO from the arguments of V\ and V 2 . 


and 


av;(* 3 ) 

dx 3 


-y/V l (x 3 )=0, x 3 =0 


dV 2 (x 3 ) 

dx 3 


+ iX~V 2 ~ 0 ’ 


x 3 ^>°° 


(23a) 

(23b) 


G(k t ,x 3 \y 3 ,co) = 


V 2 (x 3 )V 1 (y 3 )/W(y 3 ), 

Vi(x 3 )V 2 (y 3 )IW(y 3 ), 


y 3 <x 3 

y 3 >x 3 


(24) 


and W (V [ = — V,V 2 is the Wronskian. For convenience both V l and V 2 are normalized at zero 

argument y 3 = 0 such that 


Vi(0) = V 2 (0) = l, 


(25) 


and since according to Abel’s formula the Wronskian to equation (22) is independent of y 3 , then we evaluate 
W = W a at the surface. Equations (23a) and (25) conclude that 

W o (k t ,co,Z) = V 2 (0)-y/(k l ,co,Z). (26) 

The above expression shows that the Wronskian W 0 depends on the BL velocity and temperature profiles at y, = 0 , 
as well as the surface impedance Z through parameter If/ . It is noted that equation (22) supports Helmholtz 

instabilities that would arise when k, and CO satisfy the Eigen-wave condition W a (k t ,OJ,Z ) = 0 . The instability 
waves could be triggered when the BL profile supports the Eigen-wave condition. These waves could grow 
exponentially and dominate the region they occupy. Consequently, we need to set the wave number limit to avoid 
these special cases. Jones [21] studied the instability of a two-dimensional shear layer using a model velocity profile 
in which the velocity increased linearly from 0 to U e over a distance of h from the surface and then remained 
constant. He showed that the onset of Helmholtz instabilities occur at a Strouhal frequency of COh !U e ~ 1 / 3 
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where CO denotes the frequency of sound waves falling on the shear layer. At frequencies above this critical value 
the layer remained stable. 

Upon placing (24) into (10), the Fourier transformed GF at an observer point above the source region is given as 


G(k t ,x 3 \y 3 ,co) = 


i 1 -co + k l U(x 3 ) V 2 (x 3 )V l (y 3 ) 

(2k) 3 c(y 3 )c(x 3 ) (- a ) + kjJiyJf 


x 3 > y 3 


(27) 


An alternative form of the GF valid for X 3 <y 3 (see equation 24) should be used if one were interested in the 
surface pressure. The two functions V t are determined numerically. The first variable is evaluated when equation 
(22) is solved as an initial value problem with V^O) = 1 and V/(0) = Iff . The second variable is determined when 
(22) is solved as a boundary value problem subject to V 2 (0) = 1, and V 2 (x 3 ) + Z^V 2 ( x 3 ) = 0 as X 3 — > 00 . 
Fortunately, as we shall see next, there are only certain values of wave numbers k x and k 2 that make the major 
contribution to the radiated far-field noise. 


3. Stationary Phase Solution 

We apply an inverse FT with respect to wave numbers k x and k 2 and write the GF as 

G(3c,y;fi))= J J G(k t ,x 3 1 >- 3 ,C0)e nk < ( A| “ V| ( ^ >l dk t dk 2 . 


(28) 


Since the interest is in the far-field noise, the appropriate form for V 2 is 
V 2 {k t ,x 3 ,co)-b 2 e~ lXooX \ . 


(29) 


Complex number b 2 (k t ,CO) has a dimension of unity, and should be a known parameter once the boundary value 
problem (22) is solved for V 2 . Placing (27) and (29) into (28), the double integral involving wave numbers k\ and 
k 2 becomes 


G(x,y;«) = 




-(0 + k x U(x 3 ) b 2 V t (k r »y 3 .o>) 


^—e 


iQ(k t ,x,co ) 


(2k) c(y 3 )c(x 3 ) J k J{-(o + k l U(y 3 )) W 0 (k t ,co,Z ) 


dk\dk 2 


(30) 


where 


Q(k,,x,co) = k l (x 1 -y l )+k 2 (x 2 -y 2 )- x^x 3 , (31) 

and was defined in equation (12b). We adopt a spherical coordinate system where 6 is the angle the observer 
position vector makes with the direction of the mean flow (Fig. 4) 

(x l -y l ,x 2 -y 2 ,x 3 ) = R(sm(pcosd, cos0, sin^sin#), (32) 
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Figure 4. Spherical coordinates 0 <0 <n, 0 <(j)<n . 


where 7? = \ x — y | . The wave number vector in expressed in a similar fashion as 

(k l ,k 2 ,k 3 ) = K 0 (sin(l)'cos9', cos0', sin^'sin#') . (33) 

Using (32) and (33) in (31), the phase factor becomes 

Q(k t ,x,co) = Rk o ( sin < p cos 9 sin < } Y cos 9' + cos 0 cos < j Y - sin (j) sin 9 S(9' ,(f)' , M^)) , (34) 

where real function S is taken as the negative root of S 2 for those values of (6' that lend to a positive S 2 
S 2 (O' ,<p',MJ = (l-(l-M 2 J cos 2 O') sin 2 f - 2M„ cos 0 ' sin f . (35) 

Note that far-field radiation condition (14) combined with (33) also leads to the requirement that S 2 > 0. 

The large parameter in applying the stationary phase technique [22] is RK o 1 . This requires the source to 
observer distance R to be much larger than the wavelength of the radiated sound. If we define a Strouhal frequency 
as St — C0S o HJ e , where 8 0 is the BL thickness (distance from the wall at which the mean velocity is 0.99 U e ), then 

the large parameter requirement implies St » C oo S 0 / ( U e R ) which is not a very stringent condition even at 
relatively low frequency. 

The Jacobian in making coordinates transformation from ( k { ,k 2 ) to ( ) is K 2 q sin 0' sin 2 (/)' or 

dk { dk 2 = k 2 o sin 0' sin 2 ty'dO'dq )' . 


Using the new coordinates, the point of stationary phase, denoted as (6 s ,(/) s ) , is obtained when we set equal to zero 
the two partial derivatives of phase factor 0 




(36) 


It can be shown that the stationary-point angles are related to observer angles (0,0) as 

tan 9 = -S(d s ,<p s ,MJ/(M oo +( 1 - M 2 ) cos 9 s sinf ) , 
sin9 tan(f) = -S(9 s ,(j) s ,M^)/ cost/) 5 . 
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(37) 

(38) 


In the absence of flight effect we find 


0 s = 0, <t> s =$, (M„ = 0). 


( 39 ) 


Equation (39) simply implies that when flight effect is small, the major contribution to the radiated sound field is 
due to wave numbers with magnitude K o in the same direction as the observer point ( 0 , 0 ) . When 0 ^ 7T / 2 the 

two equations (37) and (38) need to be solved numerically using Newton-Raphson iteration method to determine the 
stationary point angles as a function of the observer angles and flight Mach number. In the special case when 

0 — 7T /2 equation (38) shows that (j) s = 7T / 2 , and from (37) we find 


COS0" = 


1 


f 




-M + 


cos 6 


V 


Vl-Mlsin 2 ^ 



(40) 


Figure 5 shows the range of polar angle 6 that are subject to acoustic radiation from the surface at flight when 
0 = K /2 . It is seen that the minimum value of 0 s as well as the maximum value of 6 are affected by flight. 
When 0 ^ K / 2 the stationary-point angles vary with both observer angles. A sample example is illustrated in Fig. 
6 at a selective value of 0 = 7T / 4 . 

The contribution to the integral in equation (30) in the vicinity of the stationary point is obtained when we expand 
phase factor 0 (to second order) about this point and evaluate the rest of the integrand at the point of stationary 

phased =(k;,k s 2 ). 


k s t =/e o (sin0*cos0\ cos0"), 


(41) 
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Figure 6. Stationary point angles 6 s and in flight at 0 = 71 / 4 . 


iQ(k t s ,x,co) • 2±s 


G (x,y\co) — i 


sin0" sin 0" b 2 (k t \co)V 1 (k t \y 3 ,co) (1- M „ sin0"cos0") 3 


(2 nfR cic(y 3 ) W 0 (k t s ,co,z) 


U(y 3 ). 


v 2 ’ 


1- ~ vy3/ sin cos 6 s 


as R — > °° , and 


3 =n 


exp 


— (Aft 2 + Azfi 2 +2A 3 aP) dadp 

V 2 J 


(42) 


(43) 


Parameters A. (/ = 1,2,3) are coefficients in the second-order expansion of the phase factor about the stationary 
point 


A 


1 d 2 @(k t \x,a}) 
Rk o dO s2 


1 3 2 0 = 1 3 2 P 

Rk 0 d<p s2 ’ 3 ~ Rk o d6 s d(p s 


(44) 


and Q(k t s ,X,C0) is evaluated from (34) and (35) subject to the condition stated earlier for S. General expressions 
for coefficients A. are provided in Appendix C. It is seen that the coupling term A 3 vanishes when (j) = 71 / 2 . 
The double integral in (43) is evaluated in closed form as 


3 = 


2 K 

i =exp 

^1 A,A 2 - A 2 1 


/ a A 2 \ t a \ 

- — sgn(A 1 — -f-) - — sgn(A 2 ) 
l 4 A 2 4 ) 


(45) 


In the absence of flight it is readily shown that 3 = —2 Ki / sin 0 , and 0 = k R . Subsequently 
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Jk.R 


G(x,y;«)~- 


sindsiruj) b 2 (k t ' ,(0)V t (k t s ,y 3 ,(0) 


1 


(2k) 1 R clc(y 3 ) W 0 (k t s ,co,z) 


1 - 




A 2 


, (M co = 0) (46) 


sin (j) cos 6 


where = /C o (sin0COS0, COS0) . Note that when /c, = K 0 sin0COS0 , the upper limit placed on /c, in section 
2 will be satisfied at all subsonic conditions U ■ : /c oo < 1 . At supersonic conditions only certain observer locations 
will be subject to the critical layer singularity. 

The GF G(x,t\y, T) is now obtained as 


G(x,t \y,r) = J G(x,y;co)e- ,0,< '- T) d(O . 


(47) 


Appendix D shows that equation (42) is also applicable to the acoustics of an unbounded medium (in the absence of 
a scrubbing surface) provided that the two linearly independent solutions V, and V 2 to equation (22) are evaluated 

subject to appropriate boundary conditions at X 3 — > +°° . 


4. Far-field Acoustics 

The spectral density of the far-field sound is a FT of the pressure auto-covariance 


°° r i t 

p 2 (x,co) = (Yp ) 2 j e ,m — | K\x,t)n'(x,t + T)dt 


dz . 


where T denotes some large time. Upon placing (5) and (47) into (48), it is shown that [23] 


(48) 


p 2 (x,co)=jj J G*(x,y-ll 2 \( 0 )G(x,y+ll 2 \( 0 )q(y,l,z)e imT dzdldy , (49) 

y l t=-° o 

where q denotes a two-point space-time correlation between noise generating sources at points A and B separated by 
space £ and time T , and superscript * denotes a complex conjugate. Using the source term in equation (1) we have 

j T 

q(y,l,z) = (ypf — J T(y-ll2,t)T(y+ll2,t + z)dt, T -> °o . (50) 

2 T * T 

Since the exponential phase factor in the GF equation (46) depends on the distance R =1 X— y I between observer 
point X and source point y , the product of the GF and its conjugate is evaluated at the center of the correlation 

multiplied by a phase factor cxp(— ik.^ ) , where wave number k is directed as X—y and has a magnitude of K o 

p 2 (x,co)= j\G(x,y;co)\ 2 j j q(y4,z)e im ~ ik ^dzd^dy . (51) 

y l *=-°° 
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As is a common practice in acoustic analogy noise-prediction models, the space-time FT of the source correlation 
function, denoted as the inner double integral in (51), is evaluated in closed-form once an acceptable physics-based 
source model is obtained. The modeling becomes much less cumbersome if the spatial derivatives present in the 
source were to be transferred to the GF within equation (5) and prior to forming the sound spectral density function. 
Consequently, the GF could be subject to new spatial derivatives at the source point y while source F(y,^) would 
be free of similar derivatives. 


Since the GF depends on both V\ and V 2 , which in turn depend on CO , 0 , and 0 , the pair (V l ,V 2 ) is determined 

numerically per pre-selected observer angles and frequency using the local BL mean velocity and density profiles. 
The integral over the source region in equation (49) is then evaluated to obtain the far-field noise as a contribution 
from independent correlation volume elements that comprise the source. As pointed out earlier, the volume 
integration is limited to regions with relatively stronger source strength while the mean flow is considered as locally 
parallel in both x\ and x 2 directions. 

In the limit of very low Mach number (M <0.10), it is argued [12] that there is a direct correspondence between 
the acoustic domain of the wall-pressure spectrum and the far-field spectral density of the radiated sound. This 
approximation ignores mean flow refraction and sound produced by the turbulence outside the boundary layer at 
distances that are not small relative to the acoustic wavelength. Following this argument, at large distance R from a 
hard wall region of area A, Howe [24] proposed an expression for the sound spectral density in terms of the so-called 
blocked pressure, which in the present notation is 


P 2 (x,cq) = Ak 2 0 


sin 2 6 sin 2 (j) 
R 2 


P(k t ,co) , 


(52) 


where P(k t ,CO) is the wall-pressure wave number-frequency spectrum, and is defined as a FT of the space-time 
pressure correlation function q(^ lr) ^ 2 ,z) of the wall-pressure 


I 00 

P(k t ,co) = - 2nf J* j^i^O^p(i(0T-ik t 4,)dTd^d^2 • (53) 

Here £ is the spatial separation vector of the correlation on the surface and k t is the corresponding wave-number 
vector. The wall-pressure correlation function over a rectangular area A is defined as 

^ ^ T L x L 2 

q{% i^2>0=lim— -J J j p'(y 1 ,y 2 ,0,t)p'(y l + ^ l ,y 2 + ^ 2 ,0,t + z)dy 1 dy 2 dt. (54) 

^°°Zi A _ r _ Li _ L2 

Various models for P(k t , CO) have been proposed in the literature, however there is no universally accepted form 
for the wall-pressure spectrum that could directly be substituted into (52). For example low wave number 
measurements of Sevik [13] suggest 


P(k t ,(0) __ M(U e /cJ(yJU e ) 

plvldl {co8jU e r 


24 <^<240, 0.01< — <0.15. 
U. c 


where V* ~ 0.03t/ e is the friction velocity, and 8 0 denotes the BL thickness as defined earlier. 


(55) 


More complicated expressions have been suggested for the wall-pressure wave number-frequency spectrum 
P{k t ,C0) with additional parameters that are meant to account for the attenuation and refraction of sound in the BL 
[25] and/or surface roughness [26] and dependence on the individual wave number components on the surface. The 
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validity of these empirical models remains to be confirmed experimentally. A comparison of semi-empirical 
equation (52) with the analytical solution given in equation (51) shows that the latter should account for all of the 
above-mentioned effects when source information as well as the GF are available within and outside the BL. Ideally, 
a RANS solution to the wall-bounded flow would provide both the mean flow and turbulence information, such as 
turbulent kinetic energy and its dissipation rate, required in evaluating the GF and modeling the source 

autocovariance q(y ,T) in Eq. (51) within the boundary layer as well as the main jet. 


5. Two-Dimensional Approximation 

When the two-dimensional flow is infinitely long in transverse direction x 2 (unlike in section 3 where the 
rectangular flow was limited to a finite x 2 domain), the analysis is carried out in two dimensions. Following the 
steps outlined earlier, the new GF is obtained from an equation similar to (30) 


G(x,y;«) 


(2nf c{y 3 )c{x 3 ){(-( 0 + k ] U{y 3 )) 2 W 0 {k„( 0 ,Z) 


(56) 


where phase parameter is now defined as 

0(&i ,x,CQ) = k 1 (x l -y l )~ x„x 3 , xl = (~ K o + ■ 


(57) 


It is readily shown that the point of stationary phase is k[ = K o COS 6 s , where stationary angle 6 s is solved from 
Eq. (40) as seen in figure 5. The GF is then written as 


where 


G (x,y;co) — i 


e m4x^ sin fl* b 2 (kj wW^kj ,y 3 ,co) ( l-M^cosA*^ 
(2 nf^Jclc(y,) W\k;.o,,Z) Q.WyAaaff/ 


3.-J 


exp 


-A x a} 


da - 



-in I A 


(58a) 


(58b) 


Coefficient is evaluated as before (see Appendix C), but with (f) = (f) s = K / 2 . In the special case when 
M ^ =0 it is seen that = 1 , ,X,CQ) = K q R , and k[ = K Q COS 6 . Accordingly 


G — 


(Inf 2 dc(y,) W„(k[, W.Z.) (1 _«Z3l cos8) a ’ 


(M.= 0). 


(59) 


The source correlation now uses £ = (^ 19 ^ 3 ) as a separation vector, and the far-field sound is evaluated using 
equation (49) with the volume element as dy — dy { dy 3 , and per unit flow thickness in x 2 direction. Note that both 
3D and 2D expressions for the GF could be written in a similar from if equations (46) and (59) were normalized 
with respect to their corresponding free- space GF. 
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6 . Numerical Results 

Sample Green’s function calculations are presented using analytical representations of the jet mean velocity 
and temperature profiles in the proximity of a flat surface. 

The mean velocity profile adopted here depends on the normal direction y 3 only, and is modeled for a jet with 
thickness D . = 2 .0 inches 


U(rp 

u j 


Dr\ 

tanh(^— ), 


1,, t/^ 1 „ U„ s , 1 

— ( 1 h ) +— (1 ) tanh — 

2 Uj 2 Uj d 2 




1/2 77 -I 

77-1 1/2 


77 <1.05 


77 >1.05 


(60) 


where n — V 3 / Dj is the non-dimensional distance from the wall. At 7] > 1 .05 the above mean velocity decay is 
compatible with the well-known similarity rules of Tennekes and Lumley in the parallel flow approximation [27]. 
Using parameters d 2 ) = (0.10, 2.0) , and in the absence of flight, the velocity profile divided by jet velocity 

U j is plotted (Fig. 7). The boundary layer thickness for this flow is 8 0 I D .—\ .324 d x . It is noted that at the 

connection point (7] = 1.05), the left-hand derivative U'(j /) is zero, but the right-hand derivative is slightly 
different from zero, accordingly an interpolation of second-order or higher may be required in the vicinity unless 
this point is bypassed in the integration. 

The mean static temperature is modeled as a composite of the two profiles - the first profile 7^(7]) is obtained by 
placing the above velocity profile in Crocco-Busemann law 


TM 

71 


= 1 + (?W) 


u (77) 


y-lT U(77)Y 

2 l c ~ V 


(61a) 


where T R is the jet plenum stagnation temperature ratio. The second profile is intended to simulate frictional heat 
generation in the proximity of the wall 


t 2 0 7) 1 


11 ,1,1 ^ X 

— + — tanh — ( Dm 

2 2 d,D i 77 jU 


(61b) 


Using parameters (d 3 , d 4 ) = (4.0, 3.0) the temperature profile T(j]) = T l (j])+T 2 (j\) and the sound speed 
c = are shown in Fig. 7 subject to jet exit values of U - / C m = 0.90 , and T R = 3.0 . 
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Figure 7. Profiles for mean axial velocity (solid line), static temperature 

(dashed line) and sound speed (dotted line) at Uj /c 00 = 0.90,r fi = 3.0 . 


The GF is solved numerically at specific values of the Strouhal frequency St 0 defined as 

1 0)D i 


St, = 


2n U : 


(62) 


We define non-dimensional parameters K = K n D. = 2 K St g (Uj / C, ) , c=c/c m , U = U / c„, R — R / D ■ , 
and divide the 3D and 2D expressions in (46) and (59) by their respective free-space values of(— e xR / AkR) and 

(_ e iKR+i7rl4 / 

). Additionally constant \/{kc is factored out in order to define a normalized GF 
applicable to the far field. In the absence of flight effect we have 


G N (rj,co) = 


sin 0 sin 0 b 2 (k° ,(Q)V x (k° ,T], cq) 


1 


c W 0 (k, s ,co,Z ) 

sin0 b 2 (ki ,a))Vi(ki ,ri,(D) 


(l-C/sin0cos0) 
1 


d j 


W o (k;,0),Z ) (l-[/cos0) 2 


[3D] 

[2D] 


(63) 


Here W 0 = DjW 0 is dimensionless, and derivatives present in the Wronskian are now evaluated with respect to the 
normalized distance T] . Factor ( Dj I t Q ) present in the 2D solution implies that the pressure field is evaluated per 
unit flow thickness t in X 2 direction. Since the Wronskian in equation (22) is independent of normal distance to 
the wall, computationally it is advantageous if it were evaluated at infinity where V 2 may be eliminated from the 
solution 


b 2 (k t s ,a))_ 
W 0 (k t s ,co,Z ) 


e iD ^n 

iDjXX (k t s , r\ ,co ) + d V, (k t s , r\ ,co) / drj 


(64) 
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Sample Computations 

We start with a rigid boundary condition on the wall (Z — > °o) and a set of values for the observer 
location and flow condition as {(j) = K /2, 6 = K /4, S7 0 = 0.25, C/y /c^ 0.90, 7^ = 3.0) . Two linearly 
independent solutions V, and V 2 to equation (22) are shown in figures 8 and 9. Solution V, is real, and possesses a 
relatively larger amplitude compared to the second solution V 2 . It is noteworthy that in view of equation (25), the 
amplitudes of V x and V 2 are entirely arbitrary, and the GF should be independent of this selection. At large 
distance 7] from the boundary, both VJ and V 2 exhibit an oscillatory behavior with a regular wavelength of 
2nl\Xoo I. 


Parameter b 2 (k f s ,CO) is complex, and is determined from Eq. (29) as 7] — > °°. Figure 10 shows that factor 
V 2 (Tl)exp(irjDj% 00 ) approaches constant b 2 shortly after 7] = 4 . The real and imaginary components of 

G n (TJ,(D) are presented in Fig. 11. It is seen that the zero intersects of G N follow those of . It is noted 

that the GF represents the combined effect of direct radiation from a harmonic point source at 7] plus reflection 

from the surface (and the adjacent flow), therefore zero intersect-points correspond to source locations where the 
two effects cancel each other out. Since the far-field pressure is evaluated according to equation (49) as the 
convolution of the source correlation function and the GF, the importance of the GF appears primarily within the 
region of nonzero sources. This, most likely, is limited to the initial several diameters normal to the wall where 
jet/wall turbulence and/or heat related sources are present. Figure 12 displays an expanded view of the GF within the 
first 2 diameters normal to the wall. The amplitude is relatively larger at Tj = 1 compared to regions within the BL, 
and the phase is constant at 1.16 Rad. 


VM 



Figure 8. Solution V' (7]) subject to the mean velocity and temperature profiles of Fig. 7 at 

(0 = n / 2 , 6 = n / 4, St o = 0.25, U. / = 0.90,7; = 3.0). 
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V 2 (tj) 



Figure 9. Real (solid line) and imaginary (dashed line) components of solution V 2 (T]) subject to the 
mean flow profiles of Fig. 7 at (0 = k/2, 6 = k/ 4, St o =0. 25, [/. Ic^ — 0.90,7^ = 3.0). 



Figure 10. Real (solid line) and imaginary (dashed line) components of constant 
parameter b 2 as rj — » oo . 


Gn 



Figure 11. Real (solid line) and imaginary (dashed line) components of the Green’s function 

G n (t],qj ) at ((/) = k/ 2,6 = K/4,St o =0.25,Uj /c^=0.90,T R — 3.0). 
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Figure 12. Magnitude (solid line) and phase (dashed line) of the Green’s function G n (T],co) 
at (0 = 71/ 2 , d = tt/4, St= 0 . 25 , t/. /<;„.= 0 . 90 , = 3 . 0 ). 

Next, the frequency is increased to = 0.50 and 1.0, while the remaining parameters are kept the same as before. 

Noticeable reduction in the magnitude of the GF (see Fig. 13) indicates that sound generation within the region of 
nonzero sources should decrease at higher source frequency provided that the source density remains constant. In the 
special case of a uniform mean flow (U — const , C — const ) , it is readily shown that the two linearly independent 

solutions to Eq. (22) are V, = COS(x 3 ^) and V 2 = exp (ix 3 %), which results in a Wronskian of W Q = i%. 

Parameter^ is evaluated according to Eq. (11) with k t = k f s , and is independent of the normal distance from the 
surface. Using that in Eq. (42) shows that the GF is inversely proportional to frequency. 


IGJ 



Figure 13. Dependence of the GF magnitude on source frequency 

(0 = 71 / 2 , e = 71 / 4 , U. !c_= 0 . 90 , T r = 3 . 0 ) . 


In the next example, temperature effect is examined at T R = 1.0 and 2.0, while the remaining parameters are kept the 
same as in Fig. 12. The temperature profile for the unheated case is shown in Fig. 14. As shown in Fig. 15, the GF 
magnitude near the wall (7] <1.1) is increased considerably (nearly four times) as the jet temperature ratio is 
reduced from 3.0 to 1.0. It is noted that according to Eq. (63), that temperature effect enters the computations 
through factors (c^ / c) as well as solution . 
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Figure 14. Mean flow profiles for mean axial velocity (solid line), static temperature 
(dashed line) and sound speed (dotted line) at Uj/c„=0.90,T r = 1.0. 
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Figure 15. Dependence of the GF magnitude on temperature 

(0 = nil, 6 = k! 4, St o =0. 25, Uj lc_= 0.90). 


The effect of the observer angle on the GF is examined at 3 polar angle of 6 = 45 °, 90 °, and 120 ° while the 
remaining parameters stay the same as before. Figure 16 shows a substantial increase in the magnitude of the GF 
within the BL ( Tj < 0.13 ) along the sideline (at 90°) as well as at up-stream angle of 120° compared to 45°. 
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Figure 16. Magnitude of the Green’s function G N at indicated polar angles at 

(0 = ^/2, St o =025 9 Uj/c„=0.90 9 T R = 3.0). 


Wall conditions may be examined by selecting numerical values for the normalized impedance function Z . The 
real and imaginary components of the surface impedance depend not only on the wall characteristics (such as 
resistance and reactance of a liner), but also on the external flow characteristics such as Mach number as well as 

sound frequency and amplitude [28, 29]. Aside from the solid boundary conditions (Z — > oo) examined earlier, 
here we choose a pair of values for the specific resistance and reactance as Z = (3.0,0) and Z = (1.0, 0.5). 
Computational results are shown in figure 17 at unheated condition T R = 1 .0 . At an observer angle of Q —K 1 4 , 

the amplitude of the GF is rather unaffected between Z = oo, and 3.0, and is slightly reduced at Z = (1 .0, 0.5) . 

The phase (shown in radians) is constant for the rigid wall, but becomes a function of the source location with a non- 
rigid type boundary condition. 


G i 



Figure 17. The effect of wall impedance on the GF - rigid wall: (solid line); Z = (3. 0,0.0) : (dashed 
line); Z = (1 .0, 0.5) : (dash-dot); (0 = tt/2, d = n!4 , St = 025 , 11 . !c„ = 0.90,7* = 1 .0) . 


22 


Flight effects are examined for a range of polar angles G at a constant azimuthal angle of (f) — K /2 using selective 
ambient Mach numbers of 0.0, 0.20 and 0.35. The 2D and 3D solutions scale differently with respect to flight 
factor. Using equations (42) and (58a), it is readily shown that the two GF solutions scale as (1 — sin0) 1 and 
(1 — sin0) 3/4 at (j) — K /2 , respectively. Here we choose the 3D solution to investigate the flight effects. The 
remaining parameters are selected as St Q = 0 .25 ,Uj/c oo =0 .90 ,T R = 1 .0 , and Z — > °o . Shown in figure 1 8 is 
the GF at the ambient Mach numbers of M ^ =0.0 and 0.35 as evaluated following Eq. (42). It is seen that the 

peak directivity angle rotates from 6 — 51° under static condition to 0 — 35° in flight at M oo = 0.35 . Aside 
from source location, directivity pattern is a function of frequency and jet temperature. 


e 
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Figure 18. GF directivity for indicated source ( 7] = 0.20,0.50,1 .0 ) at M oo = 0.0 (top); and 
M oo = 0.35 (bottom), with ( St Q = 0.25, U j / cj= 0.90,7^= 1 .0, 0 = 7T/2 ). 


Finally, we examine the above GF in a high-aspect ratio rectangular jet in isolation (in absence of a surface). Since 
the GF accounts for both the wall reflection and direct radiation from the source, we should expect a reduction in the 
magnitude of Gn in absence of a reflecting wall - although the shear layer itself could refract the sound in certain 
directions. The mean velocity selected for this exercise (see Fig. 19) exhibits a decay profile similar to that presented 
earlier in figure 14, and the static temperature is evaluated according to Crocco-Busemann law (61a) at a stagnation 
temperature ratio of T R = 1 .0 . The corresponding Green’ s function is readily evaluated as discussed in Appendix 
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D. A sample result is shown in figure 20 relative to a far-field observer positioned on the positive side (+X 3 ), and 
at a polar angle of 45°. Significant noise shielding is achieved when the source is located outside the shear layer on 
the opposite side of the observer (such as TJ = — 3 ). When the source is outside the shear layer but on the same side 
as the observer, the shielding represents a periodic pattern (standing-wave) due to direct radiation plus refraction 
from the shear layer. Within the jet plume or its shear layer, sources more distant from the observer are subject to 
more shielding, i.e., smaller GF. 

The polar directivity of the GF in an isolated jet is further examined when a point source at St Q — 0.25 is placed on 
the centerline ( Tj = 0 ). This is compared (see figure 21) with the earlier directivity predictions in the presence of a 
rigid wall at source locations of (7] = 0.20 and 0.50 ). An increase of 5-6 dB is predicted with surface present. 
In the absence of a source model, we qualitatively compare this last result with measurements [4] in an AR=8 
rectangular jet, with thickness Dj = 0.6696-in normal to surface, plate length X TE =12.0-in, and standoff h =0. The 
nozzle operating condition is set point SP07 (Mach 0.97, T R = 1.0, [A lcj= 0.90 ). At / =5000-//z , 

corresponding to Strouhal frequency of /D./tA = 0.28 , jet noise directivity measurements (figure 22) project a 
similar increase in noise above the surface as that predicted by the GF calculations shown in figure 21. 



Figure 19. Mean flow profiles in the absence of a nearby surface: mean axial velocity (solid line), 
static temperature (dashed line) and sound speed (dotted line), at U jlc^ — 0.90,7^ =1.0 . 
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Figure 20. Green’s function G n (T],co) in an isolated rectangular jet with an observer at 

( x 3 —> ), and ((f) = 7T/2, 6 = n! 4, St o =0.25,U j /c oo = 0.90,7^ = 1.0). 
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Figure 21. Green’s function directivity in a high AR rectangular jet, with/without a rigid surface, at 
indicated source locations with (0 = 7T/2, St o =0. 25 , Uj !c oo = 0.90,7^ =1.0). 
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Figure 22. Measured SPL directivity [Ref. 4] at a center frequency of 5000-Hz, AR= 8 rectangular 
exhaust with ( X TE = 12.0-in., /z = 0) and jet condition SP07 ( [/. 7^ = 0.90,T R = 1 .0) . 
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7.0 Summary 

Prediction of the aircraft noise due to the interaction of the jet exhaust with a nearby solid surface is of 
increasing interest in the design of future civil transport. Concepts such as hybrid wing body aircraft, over the wing 
engine mount, or distributed propulsion are likely to contribute to jet surface interaction noise. While nearby 
surfaces could provide significant noise shielding in certain directions, they also increase noise in other directions. 
In this study, we presented a formal solution to the propagation Green’s function applicable to a rectangular jet 
exhaust in the proximity of a flat surface. The governing equation is the compressible Pridmore-Brown equation. 
Both frictional heat generation and thermal gradients in the boundary layer are potential contributors to noise 
generation and propagation. The general expression for the far-field noise (Eq. 49) resembles that in the absence of a 
nearby surface - the major difference is concealed in the propagator, which should now be evaluated subject to 
prescribed boundary conditions on the surface. It was pointed out that aside from the scrubbing noise, scattered 
noise (or trailing edge noise) is also present on both side of the surface. A comprehensive prediction model needs to 
account for both noise components. 

Using analytical representations of the mean flow profiles, results were presented that show a parametric study of 
the propagation Green’s function at selective values of the source frequency, jet temperature, observer angle, wall 
impedance, and flight Mach number. It was shown that the magnitude of the GF decreases with increasing source 
frequency and/or jet mean temperature. The phase remains constant for a rigid surface, but varies with source 
location subject to a non-rigid surface. As expected, presence of an ambient Mach number sways the peak 
directivity angle to smaller downstream angles. Absence of a reflecting wall reduced the impact of the GF on the 
noise level by 5 to 6-dB. 

For round jets, a conformal mapping to a rectangular strip needs to be carried out. In practice, a Reynolds-Averaged 
Navier- Stokes solution to the nozzle flow in the proximity of the surface (with an appropriate turbulence model) 
provides the mean flow and turbulence information required for modeling both momentum-flux and enthalpy-flux 
source components. The propagation equation should be solved numerically at each pair of observer angles ( 6 , 0), 
and at each frequency CO subject to the local mean velocity and density profiles. In a 2D formulation, the GF is 
evaluated according to Eq. (58) while the mean flow is considered as a superposition of parallel slices in x 1 
direction and infinitely long in X 2 direction. When the locally parallel flow approximation is extended to X 2 
direction the GF is evaluated according to Eq. (42). Here, the jet may be sub-divided into elements in the span- wise 
X 2 direction in order to account for a finite wetted domain with an eventual decay of the source density in this 
direction. 

Nomenclature 

b 2 Far-field amplitude (eq. 29) 

c Sound speed 

X Phase (eq. 11) 

Dj Jet dimension 

S 0 Boundary layer thickness 

T] Normalized distance y 3 / ZU 

G Green’s function 

G, G Transformed Green’s function 

G n Normalized Green’s function 

h' Enthalpy fluctuations 

h Enthalpy 

h' o Moving frame stagnation enthalpy ( h' + V- v' / 2 ) 

Y Specific heats ratio 

F Source 
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k 

K 

K o 

M 

Wave number 

Wave number vector (k { ,k 2 ) 
Wave number magnitude ( CO lc ^ ) 
Mach number (U / c) 

CO 

n' 

P 

P(k t ,co) 

q(y,H,T) 

p 

R 

St. 

Radian frequency (2 nf) 

Normalized pressure fluctuation 
Acoustic pressure 

Wall-pressure wave number-frequency spectrum 

Source correlation function 
Density 

Distance \x — y\ 

Strouhal frequency f D . IU . 

1 

V 

t 

T 

1 R 

T 

u 

U e 

v' 

Spatial separation vector 
Wall function (eq. 21b) 

Time 

Stagnation temperature ratio 

Time delay 

Mean axial velocity 

Peak velocity in the boundary layer 

Fluctuating velocity component 

v, 

e 

0 

W o 

X 

y 

z 

Mass-averaged velocity component 
Polar angle 
Azimuthal angle 

Wronskian 

Rectangular coordinates 
Source location 

Normalized surface impedance 

Subscripts 

O 

oo 

At the surface variable 
At ambient conditions 

Superscripts 

Time average 

A 

Fourier transformed variable 
Favre average (q = pq! p ) 
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Appendix-A Influence of the Solid Surfaces 

By dividing the flow variables into their mean and fluctuating components, the Navier-Stokes equations may be 
rewritten as a set of mean flow equations, plus a set of five linear equations for the fluctuating components [18] 

L Vfl H M =/ v (x,0, (p,v = l,2,...,5) (Ai) 


where U denotes the five dependent fluctuating variables 


U M =( m i’P'e’P') 
m i = pv', i = 1,2,3 


(A2) 


(Latin indices repeat from 1 to 3, and Greek indices repeat from 1 to 5). The first three equations in (Al) correspond 
to the linearized momentum equation, followed by the energy and mass equations in the fourth and fifth place. The 
source term appearing on the right hand side of (Al) is 


f v (x,t) = ^ e ' v , J +S v4 (y-Y) e ';^, (i J = 1,2,3), v = 


(A3) 


< =-(P<v'-pv>'), 

1 (A4) 

v' 4 =(y-l)(h' + W 2 ), v'=0. 

The first three source components relate to difference between momentum flux terms and their Favre-averaged 
values, while the fourth term is related to the enthalpy flux and its Favre-averaged value, and the fifth source term 
(mass equation) is obviously zero. 

Note that both momentum variable m i and pressure variable p' are non-linear, however, in practice we set 
m j - pv' , and 


P 


(7-1) 

2 


e 


// 

kk ’ 


implies that p' e — p in the acoustic domain where turbulence is zero. 


(A5) 


The Green’s Function (GF) to set of equations (Al) with a non-zero delta function source placed in the Gth 
equation is denoted as 

huSno (*» * 1 y^) = AA* - y)S(t - r). (A6) 

The field variables of interest are obtained from equations (Al) and (A6) using the source/GF convolution integral 
Upixj) = J J g^y(x,t I y,T)f v (y,T)dydT, (A7) 

t y 


For example, when jl — 4 the far-field acoustic pressure is evaluated as 
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(A8) 


p\x,t ) « p'(x,t) = J J g 4v (x,t I y,r)f v (y,r)dydr . 

T y 


Considering that f 5 — 0 , the latter operation requires solving (A6) four times, each time by placing a delta function 
in the one of the first four equations. This process would also solve for other GF components that are not necessary 
for evaluating p' . It is computationally advantageous if the field variables were evaluated from the set of the 
adjoint equations that are governed by 

4 ? S ( £ (x, t\y, r) = d va 8{x - y)S(t - t) . (A9) 

The adjoint operator (denoted by superscript a) is readily evaluated once we multiply equation (Al) by 
g^ a J (x,t \ y ,T) , and rearrange its right hand side to appear as 

8n(x,t I y,z)L w u^x,t) = u^(x,t)L { “l g ( v a J(x,t I v , t ) + ^- ( g ( ^ ) + F ja . (A10) 

The last two terms on the RHS of (A 10) represent the bi-linear form of U ^ and , and 

Fja (x,t \y,T) = Vj(x,t)u^x,t)g^(x,t \y,z)+ u 4 gfj + u j + c 2 g 4 J ) . (Al l) 


We solve for U G (y, T) by placing (A9) into (A10), and integrating this equation within a large four-dimensional 
space-time volume (x,t) surrounding the source 


J J g ( v a J(x,t I y,z)f v (x,t)dxdt 


= Mff (y,T)+Jj[ -(u^) + —F ja 


dxdt . 


(A 12) 


Following the divergence theorem, the volume integral on the RHS of (A 12) may be converted into a surface 
integral. When there are no solid surfaces present to interfere with the sound, this integral is evaluated on the 
surface of a large 4D sphere in (x,t) domain - resulting in a null contribution due to the vanishing of the field 
variables U ^ on such a surface. Consequently (A 12) becomes 


u^(x,t) = J J g ( v a J(y,T I x,t)f v (y,T)dydT . 

T y 


(No solid surfaces) 


(A13) 


Equations (A7) and (A 13) point to the reciprocity of the GF in the absence of solid surfaces, i.e., 

gyjiyj I x,t) = g^ixd I y, r) . (A 14) 

By setting ju = 4 we evaluate the far-field pressure. It is seen that the advantage of (A 13) over (A7) is that we now 

place a source only in the fourth equation within set (A9) to compute the required GF components for pressure. 
Equation (A 13) represents the fundamental solution to the governing equations. In the presence of solid boundaries 
in a finite domain, the general solution, when formulated in terms of the adjoint GF, is complemented with surface 
integrals as shown in (A 12). The reciprocity condition (A 14) would not necessarily be satisfied when reflecting 
boundaries are present in the acoustic medium unless surface integrals are forced to disappear by imposing 
conditions on the adjoint GF. 
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To investigate the surface integrals further, we allow volume x(t) , externally and/or internally, be bounded by 
surface y(/)with a unit outward normal Tl (Fig. Al). The first term within the bracket on the RHS of (A12) is 
expressed as the sum of two terms following the Leibniz’s rule 



T- J «^S(^.ny,r>£c- J n-V s u^(x,t\y,z)ds , (A15) 

*x(0 s(t ) 


where V s (x,t) is the velocity at an arbitrary point on the surface 
s(t) . We now integrate both sides of (A15) within the time 
interval [-T, T] for some large time T. The first terms on the RHS 
of (A 15) vanishes due to the initial condition (in a remote past 
time), and the causality condition. Inserting (A 15) into (A 12) and 
applying the divergence theorem to the last term in (A 12), we 
find (upon interchanging dummy variables) 



Figure Al. Solid surfaces within the 
acoustic medium. 


i 

u^(x,t) = jjg^(y,z\x,t)f v (y,z)dydz + jdz J n j \y.u v g^(s,z\x,t)-F j ^s,z\x,t)\ds . (A16) 


T J 


-T T) 


Integral equation (A 16) represents a formal solution to the problem. The two integrals represent contributions to the 
acoustic field due to the volume and surface sources, respectively. In particular, the acoustic pressure is solved 
when jll = 4 . Moving the partial derivatives from source to the GF and/or imposing boundary conditions on the 

surface integrals may achieve further simplifications. For example, placing f y from (A3) into (A16), and using the 
identity 

g ( v l\y,T I x,t)^-e"j = ^(g^(y,z I x,t)e")-e"^-g ( v “\y, z I x,t). 


3v. 


followed by application of the divergence theorem on the first term on the right shows that 




dg! a) 


T y 
T 


dv 


^-( 7 - 1 )^' 

dyj dy 


„(a) 
M 
J 


i 

e"(y,z)dydz+ J t/rj n^is^z I x,t)e"dsdz 


(A 17) 


+ 


j dzj Hj [ VJ u v g { ^ (s,z I x,t) — F jfl (s, z I X, o] ds . 


-T 


Further, substituting for Fj from (All) into this last expression and imposing the no-slip boundary condition on 
the surface, (V- — V .).ft . = 0 , results in the following integral equation for pressure 


i 

Pe(x,t ) = - J J r JCy.T I x,t)e"(y,t)dydz+ J dz J rijg$(s,z I x,t)e"jds 

T y -T s 

T 

- J dzjn j \_p'Xsj)g%\s,z\xp)+m j (s,z)(g\ a 4 > + c 2 g\ a 4 ) )~\ds, 


(A 18) 


where 
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(A 19) 




v „(«) 




dy, 


Appendix-B Vanishing of the Surface Integrals 

Here we show that the boundary condition (21) eliminates the surface integrals in the Green’s function formulation 
of the acoustic field for bounded media. We start with equation (1), and apply Fourier transform with respect to 
variables x u x 2 and t as defined in Eq. (7) with the over hat carrot denoting a transformed variable 


d 2 $ 

dx: 


+ f(k t ,x 3 ,co)& = A, 


(Bl) 


where function /was defined in Eq. (11) and 


i}(k,x 3 ,co ) = 7r'(k t ,x 3 ,co)- 


c(x 3 ) 


-CO + k | U (x 3 ) 


, A (k,x 3 ,co) = 


iT 


c(x 3 )[-co + kJJ (x 3 )) 2 


(B2) 


Multiply Eq. (9) by 1? , and equation (Bl) by G , subtract the two expressions, and integrate the result with respect 


to X, 


or 




dx; 


dxl j 


'e**-** 6 ' 


dx , dx 


3 J 


dx 3 = J^GA-t?5(x 3 -y 3 )^jdx 3 . 


= $(k,y 3 ,co)-jG(k t ,x 3 \y 3 , 0 ))A(k,x 3 ,co)dx 3 


(B3) 


(B4) 


Upon recognizing that / dx 3 = — and dG / dx 3 = —i%JG at infinity, the upper limit of the expression 

on the left hand side of (B4) vanishes. On the right, since the GF is self-adjoint we switch X 3 and y 3 in the 

argument of G , and subsequently we substitute expressions for A from (B2), and G from (10) into (B4) to show 
that in view of (19), the right hand side of this equation also vanishes. 


Consequently we arrive at 

= *3=0. 

ax 3 ax 3 

Using surface boundary condition (18) in (B5) shows that 


(B5) 


5G 

y/G = 0 , x 3 = 0 

ox , 


(B6) 


where function Xj/^.CO ,Z) is given in (21b). 


Appendix- C Flight Effect Parameters 

Coefficients A t in equation (44) depend on the observer angles (0,0) , stationary point angles ( 6 s ,0 5 ) , and flight 
Mach number M ^ as 


31 


A 1 = sin 0 cos 0 sin p cos 6 s + sin 0 sin 6 


d 2 S 

dd s2 


02 

A 2 - sintpcosOsinp cosO s + cos<pcosp + sin0sin0^— 


(Cl) 


A 3 = sin<pcos6cosp sin0 s + sin0sin0 


a 2 s 


where 


d 2 S 

*s2 


dO s2 S(G s ,p,MJ 


d 2 S 


1 ( ,ds 

de 


de s dp 

-(t— ) 2 +(1 — M 2 )cos26T sin 2 p + M^cosd 5 sin^ 


A 

f 

J 


1 


df 2 S(6\P,MJ { ^ 


d 2 S 


dS 

-( ) 2 + (sin 2 6 s + M 2 cos 2 6 s ) cos 2p + M ^ cos 0* sin (, 


(C2) 


1 


1 

30 s dp ~ S(d s ,p,Mj{ dO s dp + 2 


+ — (1 - M 2 ) sin 20 s sin 2 p + M „ sin 0 s cos ft 


and 


dS cos 0 . .. . Jt dS cos 0 , , cos (b . , , 

= — ; sin0 sin0 , — — = — cos0 cos (p — : — sirup . (C3) 


30 s sin0 r ’ dp sin0 ' sin0sin0 

Parameter S(0 S ,p ,M x ) is evaluated as the negative root of S 2 in equation (35) when S 2 > 0 


Appendix-D Green’s Function in an Unbounded Medium 

In the absence of a boundary surface in the finite domain, the GF expression (42) is still formally applicable to a far- 
field observer once we recognize that the radiation condition requires 


Xoo — +yj Xlo as x 3 ~ ^ given that xt > 0 . The two boundary 

points ±°° are placed at the two opposite sides of the shear layer, far 
enough where the mean flow gradients may practically be considered as 
zero (Fig. Dl). As before, the two linearly independent solutions to 
equation (22) amount to an initial value problem for V ] and a boundary 

value problem for V 2 . Using notation Xoo = + \[xt> > we have 


Vj(x 3 ) = 1, x 3 ^-oo ? j = 1,2 


dv. 


-°°, j = 1 

+°°, j = 2 


(Dl) 


For an observer at X 3 — > +°° , the GF is evaluated per equation (42), while 



b 2 is defined according to equation (29). 


Figure Dl. Transversely sheared 
mean flow in unbounded media. 


When the observer is placed at X 3 — > — 00 , the GF should be evaluated for 

y 3 > X 3 (see equation 24); subsequently variable V 2 (k t s ,y 3 ,CO) is used in place of V x (k* ,y 3 ,CO) in (42), and b 2 is 
defined as 


b 2 = V l (k, ,x 3 ,co)e ix ~ x \ x 3 — ^ — 00 (D2) 
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